Mesoscopic phase statistics of diffuse ultrasound in dynamic matter 
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Temporal fluctuations in the phase of waves transmitted through a dynamic, strongly scattering, 
mesoscopic sample are investigated using ultrasonic waves, and compared with theoretical predic- 
tions based on circular Gaussian statistics. The fundamental role of phase in Diffusing Acoustic 
Wave Spectroscopy is revealed, and phase statistics are also shown to provide a sensitive and accu- 
rate way to probe scatterer motions at both short and long time scales. 
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For all waves, phase is irrefutably the most fundamen- 
tal property. On macroscopic scales, however, phase is 
randomized by multiple scattering or obscured by deco- 
herence. It is now generally accepted that a mesoscopic 
regime exists where wave phenomena persist on even hy- 
drodynamic scales. Mesoscopic fluctuations can some- 
times be long-range and non-Gaussian [H,[2|. The univer- 
sal conductance fluctuations are best known, originally 
discovered for electrons and later also observed with 
visible light 0] and microwaves [B| . In the optics of soft 
condensed matter, the existence of dynamic mesoscopic 
fluctuations has led to a new technique called diffusing 
wave spectroscopy (DWS) @. In the acoustic counter- 
part, diffusing acoustic wave spectroscopy (DAWS) 
the fluctuations of the scattered wave field are measured 
directly to probe the dynamics of disordered media. In 
seismology, the closely related technique of Coda Wave 
Interferometry Q is extending the range of applications 
being studied. 

For acoustic, seismic or radio waves, the phase can 
be easily extracted. While many applications, including 
interferometric techniques such as InSAR 0], make use 
of phase for precise measurements, the phase of multi- 
ply scattered waves has often been neglected, since it is 
generally more challenging to extract useful information 
from phase in multiple-scattering systems. Mesoscopic 
studies have revealed the fundamental relation of phase 
to the screening of zeros of random fields , but most 
of the literature has focussed on quantities such as the 
probability distribution functions of intensity, transmis- 
sion or conductance [l[, and docs not address the phase 
directly. Recent studies with microwaves infrared 
light [l2| and Terahertz radiation have explored fre- 
quency correlations of the phase. In this Letter, we 
study time- dependent phase fluctuations of ultrasound 
in a dynamic, strongly scattering medium, and examine 
the statistics of both the wrapped and cumulative phase 
evolution. This combination of theory and experiment 
reveals a deeper insight into the mesoscopic physics of 
multiply scattered waves, and explicitly shows the rela- 
tionship between the average phase evolution of a typical 



multiple scattering path - a crucial concept in D(A)WS 
modelling @, 0| - and the measured phase evolution of the 
transmitted waves. We also find that phase statistics can 
sometimes provide a more accurate method of measuring 
the dynamics than the field autocorrelation method that 
is used in D(A)WS. In our materials, the temporal phase 
fluctuations are too complex for traditional Doppler ul- 
trasound analysis. This Letter may be viewed as a way 
of overcoming these complications. 

The ability of ultrasonic piezoelectric transducers to 
detect the wave field allows the phase of the scattered 
ultrasound to be measured directly. In our experiments, 
we used a pulsed technique, so that the phase of multiply 
scattered waves along paths spanning a narrow range of 
path lengths could be investigated. For most of the ex- 
periments, the sample was a 12.2-mm-thick fluidized bed, 
containing 1-mm-diamctcr glass spheres suspended at a 
volume fraction of 40% by an upward flowing solution of 
60% glycerol and 40% water. A miniature hydrophone 
was used to capture the field transmitted through the 
sample in a single near- field speckle spot. The input 
pulses had a central frequency of 2.25 MHz (A = 0.71 
mm), were roughly 5 periods wide, and were repeated 
every 2 ms. Since the beads were in constant motion, 
the scattered signal was different for each input pulse, 
allowing the phase to be measured as a function of the 
evolution time T of the sample and the propagation time 
t of the waves. Since the sample hardly changed dur- 
ing the propagation time, the system appeared "frozen" 
to the individual pulses. At a fixed lapse time after 
each pulse input pulse, a short segment (about 4.5 pe- 
riods) of the transmitted waveform was recorded. Using 
a simple numerical technique [l4j |. the wrapped phase 
$(t) € (— 7r : 7r] and the amplitude A(t) > in each 
segment were determined as a function of time from the 
digitized field data. This technique is equivalent to tak- 
ing a Hilbert transform to produce the complex analytic 
signal A{i) exp[i(cjf + 3?(t))], where u> is the central fre- 
quency of the pulse. To achieve good statistical accuracy, 
10 sets of 8300 consecutive pulses were recorded. 

In order to gain insight into the temporal phase flue- 
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tuations of the multiply scattered waves, we examine the 
statistics of the phase evolution and its derivatives with 
time. The wrapped phase probability distribution P($), 
which gives the probability of measuring a phase $ at 
acoustic propagation time t and evolution time T, was 
found experimentally to be constant within statistical 
error, consistent with a complex random wave field de- 
scribed by Circular Gaussian Statistics (CGS) [H||. We 



have extended the theory of the phase within CGS [11 1 
to deal with the statistics of phase evolution, which in- 
volves the change in phase, or phase shift, with time. 
The joint probability distribution of N complex acoustic 
fields recorded at evolution times T$ of the sample is, 
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is the covariancc matrix If 



is convenient to use normalized fields so that C* 



(1) 
It 
= 1. 



Then, the off-diagonal elements of Gy are equal to the 



field autocorrelation function used in DAWS [7j : C^j = 
gi(Ti — Tj). For N = 2, two wave amplitudes and one 
phase can be integrated out from Eq. ([TJ) at constant 
phase difference A$(t) = $(T + r) — $(T). If we rewrap 
the phase difference into the interval (— it : n], we get for 
the probability distribution of phase evolution 
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where F = g\ cos A$. As the scatterers move for time r, 
the acoustic fields ipT and ipT+r decorrelate. This pro- 
cess strongly affects the statistics of the temporal phase 
evolution A<f>(r), with P(A$(r)) finally approaching the 
flat distribution for large time differences (gi = 0). 

The phase dynamics can be described quantitatively 
in terms of the variance of the change in phase along one 
typical path taken by the waves, (A^ rt (r)), which we 



call simply the "path phase variance" [16(. The DAWS 
auto-correlation function Q is related directly to this 
variance according to g\ (r) m cxp[— i(A0 2 ath (T))]. From 
Eq. ([2]) we can establish that the path phase variance is 
significantly different to the wrapped phase shift variance 
(A$ 2 (r)). The latter is associated with the superposi- 
tion of the waves from all paths at the detector that, 
for the phase, implies a highly nonlinear transformation. 
Yet, a universal relation (A$ 2 ) = /((A(/> 2 ath )) is pre- 
dicted, with no parameter that depends on the details 
of the dynamics. We exploit this universality below to 
find (A0p ath (r)) directly from the wrapped phase shift 
variance. Quite surprisingly, we will see later that un- 
wrapping the phase destroys this universality. 

The path phase variance can be related to the 
particle motion according to (A</> 2 ath (T)) ~ 
^nfc 2 (Ar 2 ol (r, £*)). Here k is the wave vector, n is the av- 
erage number of scattering events and (Ar 2 Gl (r, £*)) is the 




FIG. 1: (color online) The observed probability distribution 
of the wrapped phase evolution at five different time inter- 
vals t (symbols) and the corresponding theoretical predictions 
(solid curves). The dashed curves in (a) and (b) are the small 
r, small A"l> predictions (see text). The only fitting parame- 
ter (via the dependence of P(A4>) on gi)) is the path phase 
variance (A<^p ath ) at time r, which gives tdaws = 89 ms. For 
these data, the number of scattering events n = 34 ± 2. Note 
the wide variation in vertical scales from (a) to (d). 



relative mean square displacement of two scatterers sep- 
arated by the transport mean free path i* of the sound. 
At early times we expect ballistic motion, (Ar 2 ol (r)) = 
(AV^ 2 j)t 2 , and it is convenient to write (A0 2 ath (r)) = 

§ r2 / T DAWs> wner e t DA ws = 1/ \Jnk 2 (AV^) is the char- 
acteristic time scale beyond which the particle motion 
destroys the correlation of the acoustic field. 

At short times and small A$, Eq. ([2]) simplifies 
to F(A$) = i(A^ 2 ath )/[(A0 2 ath ) + A$ 2 ]3/2, show _ 

ing directly its dependence on (A0 2 ath (T)), and hence 
(Ar 2 cl (r, £*)). This expression has the same form as 
the probability distribution of the phase derivative $' 
with evolution time, P($') = \Qj[Q + $' 2 ] 3 / 2 , where 
Q = lim T ^ o (A0 2 ath (r))/r 2 = (3t 2 aws )-\ 

Figure 1 shows our experimental data for P(A$) at 
five values of t, along with fits to Eq. ((2|). The early 
times show a narrow peak centered at A<3> = 0, which 
broadens as the particles move further from their origi- 
nal positions. As r gets larger, the probability distribu- 
tion is indeed seen to approach the flat distribution [Fig. 
1(d)]. The agreement between theory and experiment is 
excellent over the entire range of phases and times, and 
for P(A4>) spanning nearly seven orders of magnitude. 
The fits provide accurate measurements of (A0 2 ath (T)) 
and hence of the relative mean square displacement of 
the particles. Alternatively, by using the universal rela- 
tionship <A$ 2 ) = J((A^ 2 ath )) (Fig. 2(a)), (A0 2 ath (r)> 
can directly be determined from the measured variance 
(A$ 2 ) - a simpler procedure than fitting P(A$). Both 
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FIG. 3: (color online) Comparison of theory and experiment 
for the phase derivative distribution functions P($' n '), where 
n = 1, 2 or 3 denotes the n th derivative of $ with respect to 
evolution time T. From Q = 40 rad 2 /s -2 , we find tdaws = 
91 ms. 



FIG. 2: (color online) (a) The universal relationship (solid 
curve), calculated from Eq. [2j between the variance of the 
phase shift along one path (j/-axis), and the measured phase 
of the transmitted field (x-axis). Dashed and dotted lines 
apply when the phase is unwrapped, making the relation ex- 
plicitly depend on the motion of the particles, (b) The rel- 
ative mean square displacement of the particles, (Ar 2 el (r)) 
(left axis), determined from the wrapped phase via the corre- 
sponding (A^p ath ) (right axis). We compare the results from 
the wrapped phase shift distribution (Q) and variance (□) 
with traditional DAWS measurements (★). The inset shows 
the effect of amplitude noise (see text). 



methods work well so long as (A$ 2 ) is less than its upper 
limit of 7r 2 /3 when the phase difference distribution has 
become flat. In Fig. 2(b), (Ar 2 el ) measured from the 
wrapped phase fluctuations and the conventional field 
autocorrelation are compared. The agreement between 
the two methods is excellent, giving direct experimental 
confirmation of the universal relationship shown by the 
solid curve in Fig 2(a). 

In cases where the noise in the measured signals affects 
the amplitude rather than the phase (eg. gain or DC 
offset fluctuations), the phase method is more robust for 
small t. This is illustrated in the inset to Fig. 2(b), 
which shows the effect of 2% random gain fluctuations in 
the field data; this amplitude noise clearly degrades the 
measurement of (Ar 2 el ) from g±, but does not affect the 
phase measurement. 

By considering the joint probability distribution of 
N = 4 fields in Eq. ([I]) and by integrating out 
one phase and four amplitudes, we have obtained an 
analytic expression for the joint probability distribu- 
tion of the first three phase derivatives with evolu- 
tion time P($', $",$"') [HI, from which the individ- 
ual distribution functions P(3>'), P($") and P($"') 



can be computed. They depend on three parame- 
ters Q,R and S, that in turn relate to time deriva- 
tives of the field correlation function gi(r) at r = 0: 
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fits to the three distributions give the values of Q, S and 
R. (Fig. [3]). These in turn provide a sensitive probe of 
the early time behavior of the particle motion, (Ar 2 el ) in 

324x 2 - 57x 4 - 3.3x 6 



powers of x = t/t DA ws- (Ar 2 el ) 



jim" . We emphasize that, by using this method, details 
about the motion up to the 6 th power in time can be re- 
trieved, which would be impossible from the conventional 
DAWS method. Figure [3] also shows that both theoreti- 
cal and experimental distributions follow an asymptotic 
power law decay with exponents —3, —2, — | (which sug- 
gests — (1 + ^) for the n th derivative). These slopes pro- 
vide a fit-independent test for CGS. 

To investigate the evolution of the phase over longer 
times, we study the cumulative (unwrapped) phase 
i> c (r), which can be obtained by adding or subtract- 
ing 27r whenever there is a jump of ±27r in the wrapped 
phase. The cumulative phase can be defined as Q C (T) — 
J„ T $'(T) dT, and is, by construction, a continuous ran- 
dom variable that is no longer constrained to the inter- 
val (— 7r,7r]. Its ensemble-average vanishes for fields de- 
scribed by CGS. For sufficiently long time intervals, we 
expect the cumulative phase shift A$ c (r) to approach 
the normal distribution with zero mean [181 ]. Its vari- 
ance is related to the cumulative phase derivative cor- 
relation function, CV(t) = ($'(T- |r)$'(T+ ir)>, 
which in CGS has the simple analytic form CV(r) = 
\(\ngi)" \n(l — g\) [H]. Fig. 4(a) compares theory 
and experiment for CV , where predictions based on a 
simple empirical crossover model for the particle dy- 
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FIG. 4: (color online) Comparison of theory and experiment 
for (a) the phase derivative correlation function CV , (b) the 
cumulative phase shift variance (A$^(r)). 

namics are also included Q, for which (Ar r 2 cl (r)) = 
(AV T 2 el )T 2 / (1 + t 2 /t c 2 ). The best fit is obtained for 
t c = 7tdawSi showing that both tdaws and t c can be 
determined from C$'. 

The cumulative phase shift variance can be calcu- 
lated from CV since (A$ 2 (t)) = 2 f Q T dv(r - v)C&{v) 
[1^ . Recalling the expression for GV(r) reveals that 
the variance of cumulative phase evolution is determined 
by (71 and its first two derivatives. This destroys the 
universal relation with the path phase variance, but at 
the same time this increases the sensitivity to details 
in particle motion at long times (see Fig. [Ha)). At 
short times, the cumulative phase variance increases as 
a power law with a logarithmic correction, (A$ 2 (t)) = 

| t V t daws[ L5 - 1h(t/\/3 tdaws)], while at long times, 
the cumulative phase shift evolves as a ID random pro- 
cess with finite diffusion constant: ^A$ 2 (r)) — > Da* t 
(Fig. IHb)). For r c > t D aws, wc find D A<S > = |r DA ws- 
The cumulative phase variance is quite different from the 
path phase variance (A</>p ath ), which has a finite diffusion 
constant only if the particles undergo Brownian motion, 
which they do not here. By comparing measured and 
theoretical (A$ 2 (r)} in Fig. Bib), the accurate value of 
tdaws = 89 ms was deduced from an appropriate trans- 
lation along the x direction. 

We have studied the phase evolution of ultrasonic 
waves in strongly scattering, dynamic media. It is impor- 
tant to discriminate the random phase evolution along 
one scattering path, usually studied in D(A)WS, from 
the observed phase evolution in a single speckle spot. 
Our experiments are extremely well modelled by circu- 
lar Gaussian statistics. This theory accurately predicts 
the behavior of the wrapped phase difference probability 
distribution, the variance of both the wrapped and cu- 
mulative phase shifts, and the phase derivative distribu- 



tions and correlation function. The excellent agreement 
of theory and experiment has allowed us to relate the 
observed fluctuations in phase evolution to the relative 
mean square displacement of the scatterers. The phase 
statistics are sensitive probes of the particle motion on 
both short and long time scales, and can provide more 
accurate information than the more traditional field fluc- 
tuation measurements. 
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